Equilibrium shape and dislocation nucleation in strained epitaxial nanoislands 
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We study numerically the equilibrium shapes, shape transitions and dislocation nucleation of 
small strained epitaxial islands with a two-dimensional atomistic model, using simple interatomic 
pair potentials. We first map out the phase diagram for the equilibrium island shapes as a function 
of island size (up to TV = 105 atoms) and lattice misfit with the substrate and show that nanoscopic 
islands have four generic equilibrium shapes, in contrast with predictions from the continuum theory 
of elasticity. For increasing substrate-adsorbate attraction, we find islands that form on top of a finite 
wetting layer as observed in Stranski-Krastanow growth. We also investigate energy barriers and 
transition paths for transitions between different shapes of the islands and for dislocation nucleation 
in initially coherent islands. In particular, we find that dislocations nucleate spontaneously at the 
edges of the adsorbate-substrate interface above a critical size or lattice misfit. 

PACS numbers: 81.10.Aj, 68.65.-k, 68.35.Gy 



The shape and size of islands resulting from a growth 
process has been a subject of intense experimental and 
theoretical studies QEHSEIEHS ^ E E3 There 
are still uncertainties as to whether the observed shape 
and size of islands corresponds to thermodynamic equi- 
librium state of minimum free energy, or whether they 
are limited by kinetic effects. In an equilibrium theory, 
the optimal size and shape result from a delicate balance 
of energy lowering through strain relaxation and energy 
cost through extra surface energy. Earlier works on equi- 
librium shape of coherent islands have used simple prede- 
fined faceted shapes in analytical calculations based on 
continuum elasticity theory 0, @- The resulting equi- 
librium shapes of 2D islands can be classified according 
to the relative abundance of two types of facets, namely 
shallow and steep facets only. The role of possible wet- 
ting films has also been recently considered for the 3D 
system of InAs on GaAs(OOl) |^|. However, in addition 
to assuming predefined shapes, in all these studies the 
role of possible dislocations in the islands has not been 
included. 

In this work, we use an atomistic model to numerically 
study the equilibrium shape of strained islands allowing 
for both elastic and plastic strain relaxation without as- 
sumptions on predefined shapes. We adopt a 2D model 
here, but extension to more realistic 3D systems is also 
possible. The reduced dimension allows to study all pos- 
sible configurations within feasible computer time. Most 
of our results here are for relatively small, nanoscopic 
islands up to a few hundred atoms in size to exam- 
ine deviations from the continuum theory of elasticity. 
In addition, we want to investigate the role of dislo- 
cation nucleation in determining the equilibrium shape 
and size of these islands. In an earlier approach, elastic 
and plastic strain relaxation have also been treated with 



a common formalism using vertically coupled Frenkel- 
Kontorova layers of finite length Q , where only particle 
displacements parallel to the substrate are allowed. How- 
ever, such model does not provide a realistic description 
of dislocation nucleation, whereas in our study, disloca- 
tions of arbitrary type are allowed in the final equilibrium 
configuration. 

In the 2D model we use for the strained adsorbate is- 
land and substrate the atomic layers are confined to a 
plane. Interactions between atoms in the system are de- 
scribed by a generalized Lennard-Jones (LJ) pair poten- 
tial ^2 U(r), modified 0, Q| to ensure that the po- 
tential and its first derivative vanish at a predetermined 
cutoff distance r c , given by 



U(r) = V(r), r < r ; 
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Here, r is the interatomic distance, e a b the dissociation 
energy, which can be different for substrate-substrate 
(ab = ss), adsorbate-adsorbate (ab = ff) and adsorbate- 
substrate (ab = fs) interactions, and r is the equilibrium 
distance between the atoms. For n — 6 and m — 12, U(r) 
reduces to the standard 6 — 12 LJ potential with a smooth 
cutoff. For most of the calculations we have chosen the 
values n = 5 and m = 8. In contrast to the 6—12 poten- 
tial, this has a slower falloff. When combined with the 
variation of the cutoff radius r c , this choice allows us to 
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study the effect of the range of the potential. The equilib- 
rium interatomic distance rg was set to values r ss = ro, rg 
and rf s for the substrate-substrate, adsorbate-adsorbate 
and adsorbate-substrate interactions, respectively. The 
parameter rf s for the adsorbate-substrate interaction was 
simply set as the average of the film and substrate lat- 
tice constants, i.e. r{ s — (rg + r ss )/2. The lattice misfit 
/ between the adsorbate and the substrate can thus be 
defined as 

/ = (rs - r ss )/r ss . (3) 

A positive mismatch / > corresponds to compressive 
strain and negative / < to tensile strain when the 
adsorbate island is coherent with the substrate. Cal- 
culations were performed with periodic boundary con- 
ditions for the substrate in the direction parallel to the 
adsorbate-substrate interface. In the calculations, the 
two bottom layers of the five-layer substrate were held 
fixed to simulate a semi-infinite substrate while all other 
layers were free to move. We checked that increasing the 
substrate thickness did not affect the results. Typically, 
each layer of the substrate contained about 100 — 500 
atoms. 

To obtain the equilibrium shape of the island for a 
fixed total number N of atoms without assuming any 
predetermined shapes, we use a systematic search ap- 
proach. Each initial coherent configuration is described 
by a set of integer numbers, rii specifying the number 
of atoms in successive layers of the island. In terms of 
these quantities, the two types of facets, considered in 
the previous works 0, |j| correspond to n,; — rt,; + i = 1 
for steep facets and n, — m+i = 3 for shallow facets. 
The only physical restriction we impose is that the is- 
land has a reflection symmetry about a line through the 
center and overhangs are not allowed. Then, for each 
initial configuration, molecular dynamics (MD) cooling 
was run to allow the system to relax and reach a min- 
imum energy configuration. The equilibrium shape for 
a given N is identified as the relaxed configuration with 
lowest energy among all the possible configurations. In 
the present case, this leads to complete relaxation of the 
interlayer bonds in the islands, while intralayer bonds 
remain strained. 

First, we present the results where the potential pa- 
rameters were chosen to be e ss = eg = £f s — 3410.1 
eV and r = 2.5478 A corresponding to Cu substrate 
[lij l . Phase diagrams for the equilibrium shapes as a func- 
tion of total number of island atoms N and lattice misfit 
/ are shown in Fig. 1(a), for two different cutoff radii 
r c = 3.82 A and r c = 5.3 A, which we call here the SR 
and LR cases, respectively. Different phases are labelled 
by the total number of layers in the island. In agree- 
ment with previous works there is a transition from 
a single layer configuration (uniform flat wetting layer) 
to an island configuration, above a critical size or lattice 



misfit. However, unlike the results from continuum elas- 
ticity theory, the phase diagram is not symmetric with 
respect to the misfit parameter and thus the behavior for 
compressive and tensile strained layers is quite different. 
The asymmetry is more pronounced for the LR inter- 
action potential, which includes significant contribution 
from next nearest neighbor atoms. 

The predominant shapes of the islands in the phase di- 
agram are shown in Fig. 1(b). We performed a system- 
atic classification of these shapes for 1 < N < 105 and 
— 10% < / < 10%. The strained film remains flat roughly 
between -5% < / < 5% (SR) and -3% < / < 7% (LR, 
up to N — 63). The region of 2D islands appears already 
for N — 11 for the largest misfits, and the extent of the 
flat film regime in the f — N plane shrinks for increasing 
N (cf. Fig. 1(a)). The shapes A, C, and D alternate 
with TV, the shapes A and C being most common while 
B is relatively rare. 

The critical size along the transition line from single 
layer to finite height island shape in our phase diagram 
can be fitted to a power law as N c oc f~ a , with a w 3.8, 
consistent with the result from the continuum elasticity 
approach Q. Similar behavior has also been found in a 
model of vertically coupled Frenkel-Kontorova chains [j| . 
For the SR potential, a rough estimate through simple 
bond counting yields a value for the ratio of the shal- 
low to steep facet surface energies r ~ 0.5. According 
to Ref. S this would imply an equilibrium shape of only 
either the steep facet with truncated top or steep facet 
followed by shallow facet and then truncated top, corre- 
sponding to the shapes A and C in Fig. 1(b). Although 
these general features agree with the elastic model cal- 
culations, there are important differences. First, there 
is compressive-tensile asymmetry in the shapes which is 
particularly pronounced for the LR potential. In addi- 
tion, for small islands (less than about 200 atoms) there 
are additional equilibrium shapes B and D in regions of 
the phase diagram in Fig. 1(a) that were not accounted 
for in Ref. 0- More recent work based on the elasticity 
theory Q finds that the possibility of shallow facets be- 
low the steep ones, as observed in shape B, should also be 
considered. However, shape D that we find here remains 
unaccounted for by the continuum elasticity theory. 

If the number of atoms N is sufficiently large, then the 
optimal size and shape for the island depend on the sub- 
tle balance between the strain energy, the surface energy 
of the island and the interface energy between the island 
and the substrate. Depending on the parameters of the 
interatomic potentials, phases other than those shown 
in Fig. 1(a) can occur. We consider here the situation 
where £f s > e ss = eg. For small misfit, complete wetting 
occurs, as expected. However, at large enough misfit e.g. 
f = 8 % and for £f s /e S s ~ 1-5, the behavior is totally dif- 
ferent. In this case, for a system of lateral size 200 and 
for total number of atoms N = 600, the minimum en- 
ergy configuration corresponds to that of a wetting layer 
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with an island of shape A with 26 layers high on top of 
it. This result also implies the existence of an optimal 
size for a given shape of the island. Our results thus cor- 
respond to the 2D version of the 3D islands in Ref. 0. 
Thus by varying the misfit and the substrate-adsorbate 
interactions, we find configurations corresponding to the 
commonly known different modes of adsorbate growth 
(Frank- Van der Merve, Stranski-Krastanow and Volmer- 
Weber), which occur as minimum energy configurations 
in our model. Detailed results for the different regimes 
will be published elsewhere (lif . 

We have also studied minimal energy paths for transi- 
tions from a flat layer to an island of finite height. For 
a fixed number of atoms N and misfit / close to the 
transition, we consider two shapes corresponding to the 
different states across a transition line from 1 to 2 ML 
in Fig. 1. Given these as initial and final states, we use 
the Nudged Elastic Band (NEB) 

method to gener- 
ate a (locally) minimum energy transition path between 
these two different shapes. We follow the similar ap- 
proach used in the study of defect nucleation in strained 
epitaxial films, introduced recently [lif. As an initial 
guess for the transition path we use a simple linear inter- 
polation. The resulting energy profile and configurations 
along the minimum energy path are shown in Fig. 2. As 
it is clear from the figure, there is a large energy bar- 
rier separating the two distinct minima corresponding to 
the two different shapes. This is generally true for any 
two states bordering the transition line in our phase di- 
agram in Fig. 1. Thus, depending on the time scale, 
the transition may occur too slowly to be observed dur- 
ing epitaxial growth. The existence of this large energy 
barrier, obtained with our atomistic model, supports the 
conclusion from the elastic theory calculations 0, that 
the transition can be regarded as first order. 

For sufficiently large islands or misfits we find that re- 
laxing an initial configuration with MD cooling already 
generates dislocations in the lowest energy state. This 
implies that the energy barrier for dislocation nucleation 
is zero or negligible above a critical value. We find that 
the dislocations nucleate from the edges of the adsorbate- 
substrate interface for an initially dislocation-free island. 
This is in contrast with the mechanism for a flat uniform 
film [J^, where dislocations nucleate from the top layer. 
This result provides a strong support from atomistic cal- 
culations for the conclusions obtained within continuum 
elasticity theory 0, 0- To better understand the dis- 
location nucleation mechanism, we consider a region of 
phase space where the dislocation is not necessarily spon- 
taneously generated, i.e., there may exist a finite barrier 
for nucleation. The transition path for dislocation nu- 
cleation is generated using the NEB approach with the 
coherent island and the island with dislocation as the ini- 
tial and final states. As a initial guess for the transition 
path we use again a simple linear interpolation scheme 
between the coherent and incoherent states. The result- 



ing energy profile and configurations along the minimum 
energy path found are shown in Fig. 3. The sequence 
of configurations along the transition path shows that 
the two dislocations nucleate successively at the separate 
edges of the island and then propagate inwards. The en- 
ergy barrier for a fixed island size decreases with misfit. 
This is consistent with experimental results |icj . which 
show that small islands are dislocation-free but disloca- 
tions appear in the island when it reaches a critical size. 

In summary, we have studied the equilibrium shapes, 
shape transitions and strain relaxation processes of 
strained, nanoscopic epitaxial islands The equilibrium 
shapes are determined via energy minimization for dif- 
ferent shapes without any preassumptions abut the pos- 
sible shapes, and they differ considerably from the elastic 
continuum predictions. The systematic transition path 
search combined with the NEB approach allow us to de- 
termine the energy barrier and transition path for the 
shape transition and for dislocation nucleation in ini- 
tially coherent islands. In particular, we find that dis- 
locations can nucleate spontaneously at the edges of the 
adsorbate-substrate interface above a critical size or mis- 
fit. Although these calculations were performed for a 
2D atomistic model, the method can also be extended to 
more realistic 3D systems. 
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Figure 1(a): Phase diagram showing island height as a 
function of number of atoms N and the lattice misfit / for 
the SR (main figure) and LR (inset) potentials. Different 
phases are labelled by the number of layers in the island. 
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Figure 1(b): Predominant island shapes in the phase di- 
agram for nanoscopic islands. Shape D is not predicted 
by continuum elasticity theory here. 
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INTERMEDIATE CONFIGURATION NUMBER 

Figure 2: First 12 particle configurations (upper frame) 
corresponding to the labelled extrema in the energy pro- 
file (lower frame). This minimum energy path has been 
generated by NEB and it corresponds to the shape tran- 
sition boundary for an 11-atom island to a flat layer. 
The final flat configuration (not shown here) is reached 
so that the remaining two second-layer atoms in config- 
uration 12 move close to the island edges and descend to 
the substrate. 
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Figure 3: NEB minimum energy profile and particle con- 
figurations for a 411-atom strained island, where a pair 
of dislocations nucleate at the edges of the island. Closed 
paths in (c) are the Burgers circuits around the disloca- 
tion cores in the final state. See text for details. 



